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Abstract. Accelerating Chaplygin gas combined with the decelerating braneworld 
Dvali-Gabadadze-Porrati (DGP) model can produce an overall accelerated expansion of 
the order of magnitude seen. Both models have similar asymptotic properties at early 
and late cosmic times, and are characterized by a length scale. Taking the length scales 
to be proportional one obtains a combined model with three free parameters, one more 
than the ACDM model, which fits supernovae data equally well. We further constrain it 
by the CMB shift parameter, and by requiring that the model yields a longer age of the 
Universe than that of the oldest star HE 1523-0901, U = 13.4 ± 0.8 (stat) ± 1.8 (syst). 
In contrast to generalized DGP and Chaplygin gas models, this is a genuine alternative 
to the cosmological constant model because it does not reduce to it in any limit of the 
parameter space. 
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1. Introduction 

The demonstration by SNela observations that the Universe is undergoing an accelerated 
expansion has stimulated a vigorous search for models to explain this unexpected fact. 
Since the dynamics of the Universe is conventionally described by the Friedmann- 
Lemaitre equations which follow from the Einstein equation in four dimensions, all 
modifications ultimately affect the Einstein equation. The Einstein tensor encodes 
the geometry of the Universe, the stress-energy tensor encodes the energy density. 
Thus modifications to G^ imply some alternative geometry, modifications to T^ v involve 
new forms of energy densities, that have not been observed, and which therefore are 
called dark energy. 

The traditional solution to this is the cosmological constant A which can be 
interpreted either as a modification of the geometry or as a vacuum energy term in 
T M „. This fits observational data well, in fact no competing model fits data better. But 
the problems with A are well known: its infinitesimally small value cannot be calculated 
theoretically in Quantum Field Theory, and if it can be calculated in string theories, 
these can be chosen in a nearly infinite number of ways, none of which have made any 
testable predictions. 

The search for alternatives to the cosmological constant model therefore goes on. 
No modified gravity models nor dark energy models have been strikingly successful in 
explaining the cosmic acceleration, except (at best) by introducing increased complexity 
or adding further free parameters. In this situation we think it may be worthwhile to 
try to introduce two modifications at the same time if it can be done economically. 

A simple and well-studied model of modified gravity is the Dvali-Gabadadze- 
Porrati (DGP) braneworld model (Dvali & al. pQ, Deffayet & al. [2]) in which our 
four-dimensional world is an FRW brane embedded in a five-dimensional Minkowski 
bulk. The model is characterized by a cross-over scale r c such that gravity is a four- 
dimensional theory at scales a <C r c H Q where matter behaves as pressureless dust. In 
the self-accelerating DGP branch, gravity "leaks out" into the bulk when a ~ r c H , and 
at scales a ^> r c Ho the model approaches the behavior of a cosmological constant. To 
explain the accelerated expansion which is of recent date (z ~ 0.5 or a ~ 2/3), r c H Q 
must be of the order of unity. In the self-decelerating DGP branch, gravity "leaks in" 
from the bulk at scales a ~ r c Ho, in conflict with the observed dark energy acceleration. 
Note that the self- accelerating branch has a ghost, whereas the self-decelerating branch 
is ghost-free. 

A simple and well-studied model of dark energy introduces into T^ u the density 
p v and pressure p v of an ideal fluid called Chaplygin gas (Kamenshchik & al. [3], 
Bilic & al. [4J) following Chaplygin's historical work in aerodynamics [5]. Like the 
previous model, it is also characterized by a length scale below which the gas behaves 
as pressureless dust, at late times approaching the behavior of a cosmological constant. 

Both the self-accelerating DGP model and the standard Chaplygin gas model have 
problems fitting present observational data. This has motivated generalizations to 
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higher- dimensional braneworld models which have at least one parameter more than 
ACDM, yet they fit data best in the limit where they reduce to ACDM. 

Here we combine the 2-parametric self-decelerating DGP model with the likewise 
2-parametric standard Chaplygin gas model because of the similarities in their 
asymptotic properties, taking the length scales in the models to be proportional. The 
proportionality constant subsequently disappears because of a normalizing condition at 
z — 0. Thus the model has only one parameter more than the standard ACDM model. 
It is a genuine alternative to the cosmological constant model because it does not reduce 
to it in any limit of the parameter space. 

This paper is organized as follows. In Section 2 we discuss the length scales and 
parameters in the DGP and Chaplygin gas models as was first done in Roos f6] and 
developed further in the references [HI E] • In Section 3 we discuss the basic DGP model 
in flat space, the standard Chaplygin gas model, and their amalgamation. In Section 4 
we discuss data, analyses, and fits. In Section 5 we discuss a constraint on the age of 
the Universe by comparing model predictions with the age of the oldest star. In Section 
6 we turn to the dynamical quantities w e s and q, and study their redshift dependences. 
In Section 7 we discuss the results and conclude. 



2. Length scales 

On the four-dimensional brane in the DGP model, the action of gravity is proportional to 
Mpj whereas in the bulk it is proportional to the corresponding quantity in 5 dimensions, 
M|. The cross-over length is defined as in Ref. [2], 

r c = MIJ2MI . (1) 

It is customary to associate a density parameter to this, 

Q rc = (2r c # )- 2 , (2) 

such that r c Ho is a length scale (similar to a). 

The Friedmann-Lemaitre equation in the DGP model may be written [2] 



where a = (1 + z) -1 , k = 87rG/3, and p is the total cosmic fluid energy density with 
components p m for baryonic and dark matter, and p v for whatever additional dark 
energy may be present, in our case the Chaplygin gas. Clearly the standard FLRW 
cosmology is recovered in the limit r c — > oo. In the following we shall only consider 
k = flat-space geometry. The self-accelerating branch corresponds to e = +1; we shall 
in the following consider only the self-decelerating branch with e = —1. Since ordinary 
matter does not interact with Chaplygin gas, one has separate continuity equations for 
the energy densities p m and p^, respectively. In DGP geometry the continuity equations 
for ideal fluids have the same form as in FLRW geometry [2], 

p + 3H(p + p) = . (4) 
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Pressureless dust with p = then evolves as p m (a) oc a~ 3 . The free parameters in the 
DGP model are Q rc and Q m = np m /HQ. Note that there is no curvature term since 
we have assumed flatness by setting k = in equation 

The Chaplygin gas has the barotropic equation of state p v = —A/p v [3j 0], where 
A is a constant with the dimensions of energy density squared. The continuity equation 
Q is then p\ + 3H (p v — A/ p v ) = 0, which integrates to 

Pip (a) = y/A + B/a 6 , (5) 

and where B is an integration constant. Thus this model has two free parameters. 
Obviously the limiting behavior of the energy density is 



f B ( B\ 1 ^ e ( B\ 1 ^ e 

Ptp{a) oc -^3- for a < ( — j , p^(a) oc VA for a > ( — j . (6) 

In models combining DGP gravity and Chaplygin gas dark energy [H [SJ [H] there are thus 
four free parameters, Q Tc , fl m , A, and 5, one of which shall be eliminated in the next 
Section. We now choose the two length scales, r c H and (_BA4) 1//6 , to be proportional 
by a factor x, so that 

B\ 1/6 



, , xr c /f = ;= • (7) 



It is convenient to replace the parameters A and B in Eq. ^ by two new parameters, 
tt A = E~ 2 K^A and x = 2J<X c { B / A ) 1/Q - Th e dark energy density can then be written 



PA a > 



fljK-^Vl + x^^a 2 ) ~ 3 • (8) 



3. The combined model 



Let us now return to Equation ^ and solve it for the expansion history H{a). 
Substituting f2 rc from Eq. (pi) , p v {a) from Eq. pi), and f2 m = {l° m a~ 3 , it becomes 



#(a) 



ft-. + a' 3 



n A ^i + x^An rc a 2 ) 



-3 



1/2 



(9) 



Note that Q rc and fi^ do not evolve with a, just like Q\ in the the ACDM model. 
In the limit of small a this equation reduces to two terms which evolve as a -3 / 2 , 
somewhat similarly to dust with density parameter a/^™ + VLAX 3 (AVL rc )~ 3 / 2 . In the 
limit of large a, Eq. ^ describes a de Sitter acceleration with a cosmological constant 

= — V fir c + a/^c + ^A- 

A closer inspection of Eq. ([9]) reveals that it is not properly normalized at a = 1 
to H(l)/H = 1, because the right-hand-side takes different values at different points 
in the space of the parameters Q Tc , Qa, and x. This gives us a condition: at a = 1 
we require that H(l) = xH so that Eq. ^ takes the form of a 6:th order algebraic 
equation in the variable x 

1/2 



x 



fi rc + ^ + fi AV /l + x 6 (4fi r , 



(10) 
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This condition shows that x is a function x = /(fi^, ^r c , &>a)- Finding real, positive 
roots x and substituting them into Eq. ^ would normalize the equation properly. The 
only problem is that the function cannot be expressed in closed form, so one has to 
resort to numerical iterations. The average value of x is found to be x ~ 0.956; it varies 
over the interesting part of the parameter space, but only by ~ 0.002. 



4. Data, analysis, and fits 

The data we use to test this model are the same 192 SNela data as in the compilation 
used by Davis & al. [10] which is a combination of the "passed" set in Table 9 of Wood- 
Vasey & al. [IT] and the "Gold" set in Table 6 of Riess & al. [12]. 

We are sceptical about using CMB and BAO power spectra, because they have 
been derived in FRW geometry, not in five-dimensional brane geometry. The SNela 
data are, however, robust in our analysis, since the distance moduli are derived from 
light curve shapes and fluxes, that do not depend on the choice of cosmological models. 

The Davis & al. compilation [JO] lists magnitudes /ij, magnitude errors A/Xj for 
SNela at redshifts Zi, i — 1, 192. We compute model magnitudes 

/i(zi, Q m , ^r c , M A ) = 5 Log [d L (zj, Q m , Q Ic , Q A )} + 25 , (11) 

where the luminosity distance in Mpc at redshift Zj is 

d L ( Zl ,n m ,n rc ,n A )= { —-^ — — , (12) 

H J xH(z) 
where H(z) is given by Eq. Q, and H = 72km/(s Mpc). 

We then search in the parameter space for a minimum of the x 2 sum 

2 _ f (M ~ ^{ z ii ^m, &<r c , &>a) \ 

x " h v A ^ ) ' 

As is well known in the ACDM model, the supernova data alone do not determine 
neither fl m nor well because they are strongly correlated. What the supernova data 
determine well is Q\ — Q m , but they have essentially no information on Q A + Q m . 

The situation here is similar: all the three parameters are strongly degenerate, and 
what is determined best is Qa — ^r c - Since no errors can be obtained because of the 
correlations, some further constraint is needed to break the degeneracy. One way to do 



0.24 - n c r 



that is to include as a weak CMB prior on fi^ an additional term in the x' 2 sum (13), 

2 



, . (14) 
0.09 ' 

This then permits to obtain error contours, and reduces the correlation coefficients. The 
value 0.24 comes from Table 2 of Tegmark & al. [13J, who obtained Q° m = 0.2391^ 
in a multi-parameter fit to WMAP and SDSS LRG data. To weaken the effect of this 
prior we blow the error up by a factor of 5. 

All calculations are done with the classical CERN program MINUIT (James 
& Roos[14]) which delivers x\ es ti parameter errors, error contours and parameter 
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correlations. We do not marginalize, but quote the full, simultaneous confidence region: 
a la error contour in the 3-parametric space then corresponds to xlest + 3.54 around 
the best value xlest 




Figure 1. The 1 a confidence region in the f2 rc )-plane from fits to SNela data 



and (i) the constraint (14) (the banana-shaped contour), (ii) the constraint (16)(thc 
elliptic contour) 



With the approximation x = 1 the best fit parameter values are 

Q° m = 0.26 ± O.f 6, Q rc = 0.82^, Q A = 2.2lt°j° , (15) 
3 d.f. (x 2 /d-f- = 1.029), exactly the goodness-of-fit of the 



with x 2 = 195.5 for 193 
ACDM model. 

In Fig. 1 we plot the corresponding la confidence region in the (O^, f2 r J-plane 
for a x 2 sum including Eqs. (13) and (14), as a banana-shaped contour. Clearly, the 
weak prior (14) has not done much to remove the degeneracy. The pair of parameters 
(0^,0^) is even more degenerate (not shown here). 

Instead of the rather arbitrary prior (14), we include a constraint from the CMB 
shift parameter R, which should not depend crucially on that it has not been derived 
in five-dimensional brane geometry. R is defined by 

* 1089 dz 



(16) 



xH(z) ' 

for which the value R = 1.70 ±0.03 has been measured [T4] . To permit comparison with 
the banana-shaped contour, we plot the SNela fit together with the shift parameter 
(and with x=l), a nearly elliptically shaped contour, in Fig. 1. 
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A. 



Figure 2. The contour delimits the la confidence region in the (f2 m , f2 rc )-plane in 
a fit to SNela data and the CMB shift parameter. In the blue region the Universe is 
younger than 12 Gyr. 



Fitting next the SNela data and R in the same manner as above, but with the 
correct x = 0.956, we find 

tt° m = 0.40 ± 0.05, tt rc = 1.1 ± 0.2, tt A = 2.5 ± 0.3 , (17) 



with x 2 = 195.1 for 193 - 3 d.f. (x 2 /d.f. = 1.027), slightly better than the fit g§ . We 
plot in Fig. 2 the corresponding confidence region. 

The improvement in \ 2 is due to the more exact value of x, mostly felt at small 
redshifts. The exact position of the ellipse is dependent on the value of x, as one can see 
from comparing the ellipses in Figs. 1 and 2. The shift in the position of the best fit from 



(15) to (17), is primarily due to the value of the new constraint, the shift parameter R. 



5. The age of the oldest star 

In every theory for a universe expanding with velocity H(z), the age of the universe tjj 
is given by 

" x da 

U 



HI) 



H(a) 



The WMAP collaboration [15] quotes t v = 13.73 ± 0.12 Gyr from a fit of the ACDM 
model. In the present model H(a) is given by Eq. Q with H(l) = xH . 




Figure 3. The contour delimits the la confidence region in the (f2 rc , fi^-planc in a 
fit to SNela data and the CMB shift parameter. In the blue region the Universe is 
younger than 12 Gyr. 



A model-independent limit of tjj can be obtained from the age of the oldest star, 
since tjj > t+, and this limit can be used to constrain any model for an expanding 
universe. Recently, A. Frebel & al. [16] have reported the discovery of HE 1523-0901, a 
strongly r-process-enhanced metal-poor bright giant star with detected radioactive decay 
of Th and U. For the first time, it was possible to employ several different chronometers, 
such as the U/Th, U/Ir, Th/Eu, and Th/Os ratios to measure the age of a star. From 
15 such chronometers the weighted average age of HE 1523-0901 is 13.2 Gyr. Leaving 
out the Th chronometers which have the largest systematic errors, the most useful 
value is £* = 13.4 ± 0.8 (stat) ± 1.8 {syst) Gyr. Here the systematic error is mainly due 
uncertainties in the U production ratio. 

The la statistical error, 0.8 Gyr, can be rewritten as a one-sided 68% confidence 
limit, t+ > 13.0 Gyr. The systematic error, 1.8 Gyr, cannot be handled by statistical 
methods, so we have to resort to a guess. We opt for constraining our model by t\j > 12 
Gyr. 

The effect of this constraint can be seen in Figs. 2 and 3. In Fig. 2 we over-plot the 
elliptical contour with the region tjj < 12 Gyr, where the Universe is too young, painted 
blue. The value of Qa along the contour is always the one that minimizes x 2 locally. 

In Fig. 3 we plot the la contour in the (fl rc , f^-plane with the region tjj < 12 Gyr 
painted blue. The value of Q m along the contour is always the one that minimizes x 2 
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locally. The Universe is too young if the Chaplygin gas acceleration dominates which it 
does at large values of Q^, or if the DGP deceleration is too weak which it is at small 
values of Q rc . 

The value of Q m affects tu in the same way as in the standard model: the expansion 
slows down for increasing values, and the Universe then is younger (blue). 



6. Effective dynamics 



It is of interest to study the effective dynamics of this model, as expressed by an effective 
density defined by 

p cS = p v -H/Kr c , (19) 
and an effective equation-of-state parameter 

w eS = -1 - . (20) 

Inserting p v from Eq. ^ and H from Eq. ^ into Eq. (19), one can take time derivatives 
to obtain p e ff in terms of p v and H. The algebraic expressions for p e jf and w e $(z) are 
readily calculable when parameter values are inserted, but too long to spell out here. 




Figure 4. The effective equation-of-state parameter w e g(z) as a function of redshift 
z at two points within the la confidence region. 



In Fig. 4 we show two curves for w e g (z) corresponding to selected points within the 
la confidence range. Both curves are computed at fl m = 0.344, Qa = 2.80; they differ in 
the values of Q rc : 1.07 and 0.90, respectively. At redshifts higher than zwl, dark energy 
exhibits phantom-like acceleration, u> c fr < —1, without phantom matter. In the range 
< z < 1 (depending on the parameter values) dark energy changes from phantom- 
like acceleration across the phantom divide, w e s = —1, to something like quintessence 
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matter, even to w e s > 0. At the present time one always has w e s > — 1 inside the 
la confidence range of the parameters, and in the future w e s(z) always approaches the 
cosmological constant value -1. 

In most of the parameter space w e g(z) exhibits one or two mathematical 
singularities, since p e g in Eq. (19) clearly can become temporarily negative. The first 
mathematical singularity develops from the peak and dip near z — 1 in Fig. 4, and are 
located where p e s changes sign. The second singularity, when present, develops near 
z = 2. 

Another dynamical quantity of interest is the deceleration parameter q = —1 — 
H/H 2 . For redshifts z ~ 4, q rs —1.3, just as at present. In the range < z < 4, q 
goes through a minimum of q ~ —1.8, and in the future it approaches q = —1. 



7. Discussion and conclusions 



We have studied a model combining the 2-parametric self-decelerating DGP model with 
the likewise 2-parametric standard Chaplygin gas model. The braneworld DGP model 
is an example of modified gravity which is characterized by a length scale r c Ho which 
marks the cross-over between physics occurring in our four-dimensional brane and in 
a five-dimensional bulk space. An example of dark energy is Chaplygin gas which 
has similar asymptotic properties at early and late cosmic times, and a characteristic 
length scale of its own. We take the two length scales to be proportional, and the 
proportionality constant subsequently drops out because of a normalizing condition at 
o = l. Our model then depends on only one parameter more than the ACDM model. 

The idea to combine the self-decelerating branch of DGP with some accelerating 
component has been addressed a few times before in the literature. Lue & Starkman [T7] 
and Lazkoz & al. [18J chose the cosmological constant as the accelerating component, 
mainly in order to explore models with phantom-like acceleration, w < — 1, for large 
redshifts, but which approach w = — 1 in the future. The effective acceleration is then 
increasing with time as the DGP deceleration vanishes, so that ultimately one recovers 
the standard cosmological constant model with all its conceptual problems. Since the 
accelerating component is a constant, it is not characterized by any cross-over scale, nor 
does this phantom-like acceleration ever cross the phantom divide w — — 1. Both these 
papers also discuss observational constraints and possible future signatures. 

Models which at large redshifts exhibit phantom-like acceleration, and at small 
redshifts cross the phantom divide, can also be obtained by replacing the cosmological 
constant above with a quintessence field (Chimento & al. |19j), or as here and with 
standard or generalized Chaplygin gas [7j in a decelerated DGP geometry. 

It is easy to explain the coincidence problem in the present model as well as in 
the plain DGP model: it is caused merely by the ratio of the scales of the action, the 
Planck scale Mpi on our brane and the bulk scale M5. These constants happen to have 
particular time- independent values which determine the DGP cross-over scale r c . 

We find that the effective EOS is phantom-like at large reshifts, then crosses the 
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phantom divide, so that w e s(0) > — 1 at the present time. In the future it approaches 
w = —1. 

Our model fits SNela data with the same goodness-of-fit as the the cosmological 
constant model, it also fits the shift parameter well, and over a considerable part of the 
la confidence range the age of the Universe is more than 14 Gyr, a constraint derived 
from the age of the star HE 1523-0901. In contrast to most other dark energy models, 
this model offers a genuine alternative to the cosmological constant model because it 
does not reduce to it in any limit of the parameter space. 

Our model should still be tested against other cosmological data, such as ISW 
data, CMB and BAO power spectra, all of which has to be derived in a five-dimensional 
braneworld cosmology. Such a derivation has been done recently |20j, so that these 
constraints can be included in the near future. 
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